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Abstract. Accretion disk coronae are believed to account for X-ray emission in Active Galactic Nuclei (AGNs). In this 
paper the observed emission is assumed to be due to a population of relativistic, non-thermal electrons (e.g. produced in a flare) 
injected at the top of an accretion disk magnetic loop. While electrons stream along magnetic field lines their energy distribution 
evolves in time essentially because of inverse Compton and synchrotron losses. The corresponding time dependent emission 
due, in the X-ray energy range, to the inverse Compton mechanism, has been computed. Since the typical decay time of a flare 
is shorter than the integration time for data acquisition in the X-ray domain, the resulting spectrum is derived as the temporal 
mean of the real, time-dependent, emission, as originated by a series of consecutive and identical flares. The model outcome 
is compared to both the broad band BeppoSAX X-ray data of the bright Seyfert 1 NGC 5548, and to a few general X-ray 
spectral properties of Seyfert Is as a class. The good agreement between model and observations suggests that the presently 
proposed non-thermal, non-stationary model could be a plausible explanation of AGN X-ray emission, as an alternative to 
thermal coronae models. 
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1. Introduction 

In these last years the idea of a hot corona lying above the AGN 
accretion disk has been developed in order to account for sub¬ 
stantial emission in the X-ray energy range. Indeed, the orig¬ 
inal suggestion of this scenario dates back to 1977 (Liang & 
Price [T977t . referring to the X-ray binaries’ context, and to 
1979 (Liang [THU, more specifically regarding AGN X-ray 
emission. The hot corona can, in fact, Comptonize the soft pho¬ 
tons emitted by the “cool” accretion disk, thus producing high 
energy radiation. Several scenarios have been developed for 
AGN coronae, both analytically and numerically (Poutanen & 
Svensson ll9961 Dove et al. 1 19971 Di Matteo [T998l Merloni & 
Fabian EUOTl Miller & Stone l200nt . Haardt & Maraschi <19911 
I1993> analyzed in detail and developed the idea of a radiative 
coupling between the two phases of the system (disk-Hcorona); 
subsequently, in order to better match the observations, Haardt 
et al. <19941 1 have developed a more detailed model of a non- 
uniform, blob-like corona. From the physical point of view, 
in the above papers the hot corona is assumed as the a re¬ 
gion where at least part of the accretion gravitational energy 
is released. The existence of an optically thin coronal structure, 
sandwiching the disk and characterized by a strong magnetic 
coupling with the optically thick disk, was firstly suggested by 
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Galeev et al. <1979> . The idea is that the accretion disk dif¬ 
ferential rotation, together with the natural presence of a mag¬ 
netic field, could reproduce a situation analogous to that present 
in the atmospheres of convective stars. As Heyvaerts & Priest 
( [T^ have shown, loops and arcades can form also in AGN 
accretion disks; these structures, connecting remote points of 
the disk itself, can convert disk kinetic energy into magnetic 
energy and subsequently dissipate it by emitting the observed 
spectrum. Magnetic reconnection is invoked as the physical 
process responsible for this last conversion, even though all the 
details of the process itself are not yet understood. 

An issue that has not yet been clarified is whether the anal¬ 
ogy between stellar and AGN coronae is limited to the above 
described features or whether it is more comprehensive. Once 
a physical description is chosen in which the magnetic field is 
responsible for transferring accretion energy from the disk to 
the AGN corona, we can ask ourselves whether the magnetic 
field can act in a different way and produce a different envi¬ 
ronment with respect to the stellar case. It seems reasonable 
that if the global picture of a stellar corona is applicable in the 
AGN context, the same should be true for the relevant phys¬ 
ical processes as well. The analysis of solar flares has shown 
(see Masuda et al. <1994> for a short summary) that acceler¬ 
ated particles are responsible for the hard (> 10 keV) X-ray 
emission via the bremsstrahlung interaction with thermal mat¬ 
ter. Moreover, even in the absence of conspicuous flares. X-ray 
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brightenings and type III radio bursts show the presence of non- 
thermal electron beams in the solar corona (Klein et al. ll997t . 
Hence, the existence of impulsive injection of non-thermal par¬ 
ticles seems to be a common ingredient in the physics of stellar 
coronae. The question we are interested in is whether, in the 
specific AGN context, an impulsive, non-thermal electron dis¬ 
tribution can be responsible for the emission of the observed 
X-ray radiation or, on the contrary, a stationary and thermal 
electron population at high temperature (like the one adopted 
in the model proposed by Haardt & Maraschi J1991ll993t . and 
Haardt et al. Comptonizing the soft disk photons is the 

only possible explanation of the observed high energy emis¬ 
sion. 

Stationary non-thermal models have been developed start¬ 
ing from the 1980s in order to explain AGN X-ray emission 
(for a review see S vensson i 1 994l and Ghisellini ll994t : however, 
these models, at least in their “standard” versions (Svensson 
119941 Svensson 1 1996t . showed some problems as for the re¬ 
production of high energy observations of Seyfert spectra by 
OSSE, in particular referring to the indications of high energy 
cutoffs in the spectra, and to the lack of detections at ~ 0.75- 
3 MeV by COMPTEL (Svensson [T993 Haardt fTWTl . In our 
framework, impulsive injection of the relativistic electron pop¬ 
ulation is determinant for what regards the specific character¬ 
istics of the non-thermal model we are proposing, because it 
is exactly this feature that makes the model intrinsically non¬ 
stationary. This peculiar property of our scenario distinguishes 
it from previous “standard” non-thermal models and will be 
explained and discussed in the following Sections. 

Therefore, starting from the knowledge that we have of so¬ 
lar/stellar coronae mechanisms and trying to extend it to the 
AGN coronae scenario, we shall investigate the time evolution 
(see Sect. 2) and the emission (see Sect. 3) of an ensemble of 
accelerated non-thermal electrons injected in a magnetic loop. 
How these accelerated particles are produced is not relevant to 
the following analysis and, on the other hand, it is still mat¬ 
ter of debate even in the solar corona, although it seems to be 
plausibly related to magnetic activity processes (such as recon¬ 
nection). 

In Sect. 4 we describe and discuss the application of the model 
that we have devised to the Seyfert 1 galaxy NGC 5548; in 
Sect. 5 the results of this analysis are compared to BeppoSAX 
observations. Sect. 6 is devoted to the discussion of the de¬ 
scriptive capability of the model regarding to general trends 
and properties identified for Seyfert 1 X-ray spectra, in order to 
test our model’s relevance with respect to this type of source. 
Einally, Section 7 includes both an analysis of our model as 
compared with previous models in literature, and the conclu¬ 
sions we draw from the present work. 


2. Behavior of accelerated electrons 

When a certain amount of accelerated electrons is injected 
around the top of a magnetic loop, part of them is trapped 
within the curved magnetic field lines, while part of them 
precipitates in the denser plasma, where the loop bases are 
anchored. The electrons which are magnetically reflected in¬ 


side the loop loose their energy by colliding with the thermal 
plasma present in the loop and by emitting synchrotron, in¬ 
verse Compton (IC) and bremsstrahlung (BR) radiation. How 
long it takes to thermalize these electrons and how the elec¬ 
tron energy distribution changes in time is analyzed hereafter. 
The group of precipitating electrons interacts with the denser 
plasma (of the disk or of the photosphere) and stops after a 
rapid bremsstrahlung emission. This emission is what we iden¬ 
tify with the hard X-ray (> 10 keV) emission in solar flares. 

The relative importance of the trapped and of the precipi¬ 
tated electrons depends on the specific physical configuration 
and on the density gradient of thermal matter in the loop. 

An estimate of the typical loop extension in the AGN accre¬ 
tion disk corona context can be obtained for instance following 
Di Matteo J1998> : extending the results of solar atmosphere 
magnetic buoyancy simulations (Shibata et al. I1989> . this au¬ 
thor gives an estimate of the loop top height (//flare), which 
reads //flare ~ 8//disk, where //disk is the pressure scale height of 
the accretion disk. Eurthermore, if, still following Di Matteo 
(ll998|l, we remind that, for a standard Shakura & Sunyaev 
J1973t disk, it is //disk//^ ~ 0.1 (where R is the radial distance 
from the central black hole on the disk midplane), we can relate 
//flare to the SchwarzschOd radius, Rs = 2GM^yilc^-, in fact, for 
the accreting disk central regions we can assume a represen¬ 
tative radial distance from the black hole R ~ 10/?s, and this 
leads to the following approximate relation //flare ^ 8/?s and, 
finally, to the estimate //flare ^ 2.36 x 10 '^(Mbh/10^Mo) cm. 
Other approximate estimates of coronal region extension can 
be found in the work of Miller & Stone (I2000> . whose MHD 
numerical simulations of a coronal structure formation above a 
weakly magnetized accretion disk indicate that, in the specific 
conditions of the simulation at least, most of the buoyantly ris¬ 
ing magnetic energy is dissipated between 3 and 5 disk scale 
heights (i.e., what we have called //disk) above the disk itself, 
and in the work of Liu, Mineshige & Shibata J2002t . who cite 
a representative size value ~ 10/?s cm. All the above cited esti¬ 
mates are, in conclusion, in substantial agreement with an eval¬ 
uation of //flare Tinging from ~ 3Rs to ~ IQRs- 

This allows us to finally estimate the typical electron dy¬ 
namical time in the loop, that can be defined as fdyn = //flare /c 
and represents the time it takes for the electrons to reach the 
disk; we also define a characteristic radiation time, femis, as the 
typical time it takes for the radiation emitted around 10 keV 
(the energy range we are mainly interested in) to decrease by 
a factor of twenty in luminosity. When femis < fdyn the elec¬ 
tron distribution evolves on short time scales with respect to 
the dynamical timescale and it is depleted before reaching the 
the impact region; in this limit the energy lost as impulsive 
bremsstrahlung emission can be disregarded. This appears to 
be the case in the AGN coronal environment, as it will be veri¬ 
fied in Sect. 4. 

2.1. Electron energy loss rates 

Given an electron energy distribution 

n(7o, f = 0) = noijo - 1)“^ el./cmVf/y, 
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where yo = jit = 0), in the energy range (in units of nigC^, 
where me is the electron rest mass) 

1 < Tot < To < Jo2, 


this distribution evolves subject to the system energy losses. 
The relevant energy loss rates for the present problem are in 
principle the following (see, e.g., Ginzburg 1 19691 Lang 1 19801 
Blumenthal & Gould fT^ l: 

• collision losses in the relativistic limit; 

^ = -1.5 X 10-^W[13A + In(rMh)] s-', 
at 

where nth is the thermal density inside the loop and it is ex¬ 
pressed in cm“^; this expression can be further approximated 
as follows 

— -1.5 X 10^'Vh[74.4 - In(nth)] s^* 

at 

= Ao(nth) s"', (2) 

since ln(y) is a small, slowly varying term in the limit where 
collision losses are important (see later, eq. (|^); 

• bremsstrahlung losses: 

^ = -1.4 X 10^'V/t[/n(y) -h 0.36]r s"' 
at 

= Ai(nth,y)y s^' (3) 

• inverse Compton and synchrotron losses: 

5 = -3.2 X IQ-\U + |^)y2 = A2{U, B)j^ s-' (4) 

at on 

where B is the magnetic field (in gauss) and U is the energy 
density (erg/cm^) of soft photon radiation field present in the 
loop structure. If Lsoft(f) (erg s^'keV'') is the seed (IR-UV) 
photon source spectral luminosity coming from the AGN disk 
and Liocai(e) (erg s^^keV^') is the radiation field luminosity (see 
eq. (0 in Sect. 3 for its explicit definition) emitted by the elec¬ 
trons of the relativistic distribution through synchrotron mech¬ 
anism in the local magnetic field, then U (erg cm"^) can be 
derived as 

L = ■" T'geom(^UV-x) 

47rc 

~ 17^ I ^local(^) (5) 

47rcRf„ 

where e is the photon energy, that we express in keV, in the 
range emin < e < eiviax. In this expression we take emin = 
10“^ keV and fMax = 0.1 keV, Ruv-x is the mean distance be¬ 
tween the source of soft photons and the X-ray emitting loop 
structure and Rem is the characteristic length scale of the emit¬ 
ting region. The quantity T’geom(^uv-x) is sort of a geometric 
factor, accounting for the fact that soft radiation from the disk 
reaches the coronal region from different different positions in 
the disk and therefore is differently “diluted”. We adopt a sim¬ 
plified treatment of this effect, following Ghisellini et al. (I2004t 
(see their eq. (14)), so that it can be represented by the quan¬ 
tity R'geom(Ruv-x), which depends basically only on a length 


r 


Tsoft(e) + 


scale representative of the vertical distance of the X-ray emit¬ 
ting structure from the disk itself, namely Ruv-x- 

All the above energy losses make the electron distribu¬ 
tion change its energy dependence. However, depending on the 
physical conditions of the plasma, their relative importance can 
be rather different, and in the following we discuss this issue 
referring to the coronal environment of our model scenario. 

First, bremsstrahlung losses do not seem to play an impor¬ 
tant role in the electron distribution evolution for the typical 
range of n* that we can estimate for our problem: they are neg¬ 
ligible with respect to collision losses for y < 200 and with re¬ 
spect to Inverse Compton and synchrotron losses for y > 200, 
if R > 50 gauss and n± < lO" cm“^. In the context of AGN 
coronae these limitations on the values of magnetic field and 
thermal density seem to be reasonably met and they do not 
pose very stringent conditions. In particular, this is true in the 
framework of the model we are defining for coronal emission, 
whose origin we imagine to be non-thermal. In fact, we have to 
require that a thermal component in the coronal loops has neg¬ 
ligible effects on the resulting X-ray spectrum. In other words, 
we want to disregard any possible thermal Comptonization or 
recoil effects on the spectrum itself. This implies a tighter up¬ 
per limit for the density (nth) of the thermal component of the 
coronal loops (see Sect. 4.3 for a more detailed discussion and 
estimate for the case of NGC 5548). In the light of the above 
considerations, we can therefore safely conclude that, in the 
context of AGN coronae and in particular in our model sce¬ 
nario, the range of plausible physical conditions is such that 
bremsstrahlung losses are basically ineffective with respect to 
those due to IC and synchrotron processes. 

Secondly, for what regards collision losses, we can note that 
at low energies these will make the distribution flatter, while 
inverse Compton and synchrotron losses will steepen it at high 
energies. A critical value for the electron energy can be easily 
derived comparing eq.s 0 and 0, when Aq > A 2 : 


riim = 

Inverse Compton and synchrotron losses are the most ef¬ 
ficient for y > yiinj, while for y < collision losses 
are the most important ones. Adopting a representative value 
for the coronal structure scale length ~ 5Rs (see Sect. 2), a 
rough estimate of the soft radiation energy density as U ~ 
0.0lLEdd/4ttc(5Rsf gives yi™ ~ 4 for Mbh = IO’Mq and for 
a reasonable (see Section 4.3) upper limit n± ~ lO" cm~^ for 
the thermal density. Furthermore, from eq. 0 it is clear that, 
when Ao < A 2 , IC and synchrotron losses will be dominant in 
any case, i.e. for any meaningful value of y. With the same or¬ 
der of magnitude evaluation of U that we have used above, we 
can define an estimate of the upper limit of the thermal density 
to maintain the condition Aq < A 2 fulfilled, so that collision 
losses are effectively “negligible” with respect to IC and syn¬ 
chrotron losses in the electron distribution evolution for any 
possible physical value of y. It turns out that, for this to occur, 
nth < 6.4 X 10 '° cm~^. Since for our model consistency (see be¬ 
low in the present Section and in Section 4.3 for more details) 
we have to require a much lower thermal density in the coronal 



U + B^l%n' 


(6) 
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loops, the discussion above clearly shows that in our context 
collisional losses can be disregarded as well. 


2.2. The evolution of the electron distribution: an 
analytic solution 


Melrose & Brown J1976> derived a formal method for the so¬ 
lution of the general equation for the evolution of a distribution 
function n{y, t) for particles in energy phase space. In the limit 
of negligibility of precipitating electrons (which applies to our 
problem as discussed in Sect. 2) this formal solution gets sim¬ 
plified. Taking also into account the further simplification due 
to the fact that bremsstrahlung and collisional losses are neg¬ 
ligible, the remaining relevant loss rates can be expressed as 
proportional to a power of the electron energy y with a coef¬ 
ficient A 2 that does not depend on y itself, but may depend 
on time through the locally produced (synchrotron) luminos¬ 
ity, Liocai- This condition proves quite useful, since it makes it 
possible to derive, following Melrose & Brown J1976> formal 
procedure, an analytical solution for the time evolution of the 
initial electron energy distribution n{yo, 0). 

The resulting time dependent electron distribution after a 
time t can be expressed as ; 


n{y, t) 


1 

[1 + y 6»]2 


y 

1 +y6 



(7) 


with ^ 

0 — I A2(t')dt'. 

Jo 

The range of y over which the distribution extends is also 
evolving in time and it becomes 


where 


yi(t) < y(t) < 72 ( 1 ) 


yi(0 = 


yoi 

1 - yo\d 


yiit) = 


yo2 

1 -JolO' 


3. The overall emission 

From the above expressions for the time dependent electron 
distribution, inverse Compton, synchrotron and bremsstrahlung 
emission can be computed. Obviously the spectral luminosity 
due to those radiation processes evolves in time as well, owing 
to the time evolution of the electron distribution. As a matter 
of fact, in our context, the IC component is the sole signifi¬ 
cant contribution to the resulting X-ray emission in our model. 
Indeed, relativistic bremsstrahlung emission is not explicitly 
described in this Section, since it turns out to be negligible in 
the present scenario, as we briefly discuss in the following. 

In fact, we have already shown, at the end of Sect. 2.1, that 
the effects of the bremsstrahlung energy loss rate on the global 
evolution of the relativistic electron energy spectrum are unim¬ 
portant with respect to those of the other loss mechanisms at 
work, for whatever value of the electron energy y, under the 
physical conditions expected for an AGN coronal structure like 


the one we want to model. This clearly implies that the cor¬ 
responding bremsstrahlung emission has to be negligible as 
well when compared to the IC component and this is indeed 
what we have verified from our spectral luminosity computa¬ 
tions. This is important to stress, since the resulting relativistic 
bremsstrahlung would peak at energies higher than those in the 
typical range of the first scattering IC component, so that, in 
principle, it could also happen that, although small, in the very 
high energy range the bremsstrahlung relative contribution to 
the total emission could be significant. 

The effects of the synchrotron emission of relativistic elec¬ 
trons in the ambient magnetic held have been explicitly taken 
into account, as far as the electron energy losses (and, as a 
consequence, the relativistic electron population evolution) are 
concerned (see eq. @) and for their contribution to the seed 
photon spectral luminosity that is to be reprocessed through 
inverse Compton mechanism. Synchrotron emission does not 
show up directly in the resulting spectrum since, for the whole 
of the ranges of physical parameters that we have selected as 
appropriate for the present context and explored by computing 
the resulting spectra, it turns out that 

a) due to the typical range of energies of the injected relativistic 
electron distribution (yMax ~ 10^), the synchrotron component 
covers a range of energies which does not belong to the X-ray 
domain, and, for reasonable values of the magnetic held, peaks 
at energies well below the optical and, moreover, 

b) its emission level turns out to be below that of the ob¬ 
served IR-opt.-UV spectrum as reported in literature (see next 
Section). 

However, due to the fact that synchrotron radiation is pro¬ 
duced in the same place, i.e. a coronal loop where inverse 
Compton reprocessing occurs, the geometrical factor to be ac¬ 
counted for in the computation of its energy density makes its 
contribution non negligible as far as the seed photon energy 
density is concerned. This contribution is taken explicitly into 
account by the term Liocai in eq. 0 and in the following eq. 0, 
which defines the resulting inverse Compton spectral luminos¬ 
ity as a function of time. 

Defining E (in keV) as the emitted photon energy and re¬ 
membering that e is the incident soft photon energy, the inverse 
Compton spectral luminosity, in erg s ' keV *, at a given time 
t reads 

Uc(E, t) = 7.5 X 10-15 ^ ^‘°°P^ X (8) 

Anc 


pfMa.x nj: 

Jy\ 


giy, e, t) 


TsQft(6) ^ 1 Tiocal(^; 0 

Tgeom ^2 15 

Kp 


dyde. 


where 


giy^ e, 0 = 



(1-5#:) 


for < 1 

Ay-€ 

for > 1. 

Ay^e 


In the above expression, Vjoops is the total volume of the part 
of the corona where flares take place (contributing to the emis¬ 
sion we observe) and Lsoft(e) is the illuminating IR-UV spectral 
luminosity from the accretion disk defined over the energy in¬ 
terval [emin, fMax] (sce Sect. 2.1). We assume this soft spectrum 
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as known, and we model it as a parameterized power-law with 
a high energy cut-off, as it is explicitly specified below: 

isoft(e) = < e. (9) 

In this expression UV2 is the photon index of the spectral distri¬ 
bution, Ecut is the cut-off energy (in keV), and Lsoft is expressed 
in erg s * keV *. Therefore, for each AGN source we intend to 
study, first we shall have to appropriately choose the parame¬ 
ters of our representation of Lsoft(e), so as to ht the observed 
IR-UV spectrum, as described in literature. 

The other term present in eq.(|^ , LiocaU takes into account 
any production of soft photons by the relativistic electron dis¬ 
tribution, namely synchrotron emission in the loop magnetic 
held. As outlined above, for our parameter ranges, synchrotron 
emission is too weak to be relevant in the observed spectrum, 
but its contribution as seed radiation for the IC reprocessing 
is non negligible. The synchrotron spectral luminosity of each 
single active loop is given by: 

Tiocai(G t) = 2.8 X 10^^-^^ B X (10) 

q 

P'2 ^ 

I dyn{y,t)— I ergs^'keV^', 

xjy\ Ueje^ 

where 

ec = 1.7x keV, (11) 

q is the number of active loops in the coronal structure (as¬ 
sumed constant in time) and ^^ 5/3 is the modihed Bessel func¬ 
tion of order 5/3. 

Expression 0 for the resulting IC spectral luminosity takes 
into account only the hrst inverse Compton scattering of the 
soft photons illuminating the loop, in the hypothesis that the 
second and the higher order ones have negligible effects on the 
resulting high energy spectrum in our context. The validity of 
this approximation will be discussed and verified in Sect. 5. 

It is also important to stress that, even if this model does 
not belong to the “family” of thermal models, it shares with 
them the existence of a feedback between the flaring corona 
and the cold underneath disk. In fact, we assume that roughly 
one half of the inverse Compton emitted radiation is directed 
inwards, to the disk surface; this approximate fraction of the 
X-ray radiation will therefore be somehow reprocessed by the 
disk material (absorbed or “reflected”), thus implying some de¬ 
pendence of the disk physical conditions on the coronal emis¬ 
sion and possibly some sort of relation between the high energy 
coronal emission itself and the soft disk emission. The expected 
complex interplay between disk and corona has been widely 
discussed by Haardt et al. (fT^ . However, it is important to 
stress that, in accordance with Uttley et al. conclusions, 

in this work we do not demand that soft radiation from the disk 
is only due to reprocessing of the X-ray coronal emission. 

4. Application of our model to AGN X-ray sources: 
the example of NGC 5548 

The choice of adequate observational data for testing our 
model is based on several criteria. Seyfert 1 galaxies (Sy 1) are 


excellent candidates because a substantial amount of their ob¬ 
served X-ray emission is believed to originate directly from the 
inner central regions of the accretion disk (with some repro¬ 
cessing in the form of reflection and absorption). The contribu¬ 
tion from the jet to the observed flux is believed to be negligible 
in this type of AGN. 

To exemplify the applicability of our model and to fully 
illustrate our procedure in its details, we have to choose a par¬ 
ticular source and test the model itself against the source’s 
observed X-ray spectrum. This is the subject of the present 
Section. However, the plausibility and validity of the emission 
model we present must be evaluated on the basis of a more 
general comparison with the observed properties of X-ray spec¬ 
tra of Seyfert 1 active nuclei. This comparison is postponed to 
Section 6 , in which the descriptive capabilities of our model 
will be analyzed in order to account for a few general observed 
trends and properties of Seyfert 1 X-ray spectra. 

For the present purposes of quantitative and detailed com¬ 
parison between theory and observations, it is important to 
choose a bright source for which broad-band X-ray spectra are 
available as well as good optic al/UV data. For the X-ray data 
the choice of BeppoSAX is obvious, because of its unparal¬ 
leled broad band and sensitivity (Boella et al. ll997al . Ideally, 
the optical/UV and X-ray data should be simultaneous with the 
X-ray data, however such data sets are extremely rare. 

NGC 5548 is a well-known and nearby (z=0.017) Seyfert 
1.5 galaxy, bright in the X-ray domain. It was observed several 
times by BeppoSAX (see Nicastro et al. 120001 Petrucci et al. 
120001 Pounds et al. 120031 and a vast amount of optical and 
UV observations are available for it. It is among the 8 brightest 
type 1 Seyfert galaxies observed by BeppoSAX (^ 2-10 keV ^ 
4 X 10“*' erg cm“^ s"'). 

BeppoSAX observed NGC 5548 at 3 different epochs: 
August 1997, December 1999 and July 2001. The source was 
at it brightest during the 1999 and 2001 observations. The July 
2001 data set was selected for our study because the exposure 
at low energies (LECS instrument) is better than in 1999. We 
do not have simultaneous multi-wavelength data for this source 
(extending from the IR to the X-ray range) and, as a conse¬ 
quence, we rely on literature knowledge of the IR-UV spectrum 
in what follows. 

4.1. The adopted soft photon spectrum 

The parameters describing its IR to UV spectral luminosity in 
the simplihed representation that we have introduced for this 
quantity in the previous section (see eq. ®) have been selected 
by analyzing the results described in the literature as follows. 
At energies higher than emin = 10“^ eV (124 //) we assume 
UV2 =1.3 and a normalization luminosity, Lsoft(2.25 eV) = 
2.8 X 10^® erg s“' keV'*, derived from Ward et al. (inH3. 
Since the source is highly variable this value for the luminosity 
seems a good compromise between the higher value quoted by 
Malkan & Sargent J1982I I. Lsoft(2.25 eV) = 3.6 x 10^^ erg s“* 
keV“*, and the luminosity range reported by Wamsteker et al. 
(f!^ . Lsoft(2.25 eV) = 1.3 A 2.9 x 10^'’ erg s^’ keV"'. The 
value Ho = 50 km s“'Mpc“' has been assumed to derive the 
above luminosities. Note that between emin and ei = 1.8 eV 
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Fig. 1. Dotted curves show the inverse Compton spectrum ob¬ 
tained at different times (to - 0.02 s, fi = 1.0 s, = 5.0 s 
after the injection of an electron distribution characterized by 
7oi = 50, 7 o 2 = 1000, i = 3.0 and WoFioops = 2.57 x 10^^ 
el/dy (see the text for this parameter’s definition). The contin¬ 
uous curve drawn for energies £<0.1 keV shows the adopted 
representation of the soft photon spectrum for the case of NGC 
5548. The other continuous curve, extending up to very high 
energy, is the spectral energy distribution of the time averaged 
resulting emission, as derived from eq. (O. 

(0.7 fi) the observed spectrum, like the one shown in Fig. [2 
can be represented by a power-law with photon index UV 1 = 2 
(Carleton et al. fT987t . 

An exponential cut-off at Fcut = 0.015 keV, modifies the 
slope of the UV spectrum at high energies. The value has been 
chosen so that the energy spectrum decreases for A < 1060 A 
in accordance with Brotherton et al. i2QQ2i . In this way the 
resulting seed spectrum is also in accordance with the EUV 
data reported in Chiang et al. <20001 1 where Leuv(0.163 keV)= 
8.2x10“*^ H- 9.5 X lO^"* erg s^' keV^'. 

4.2. The definition of our X-ray model spectrum: an 
illustration 

In the following, unless otherwise specihed, we translate the 
results of our calculations of spectral luminosities into the cor¬ 
responding spectral energy distribution Ef(E) (keV cm“^ s“'), 
which is commonly used and directly comparable to observa¬ 
tions, through the obvious relation Ef(E) = EL(E)l(AnD^), 
where D is the distance to the chosen source, appropriately 
evaluated. 

Fig. [H represents an illustrative summary of the results of 
our model, as it will be extensively explained in the follow¬ 
ing. In fact, this figure shows different curves referring to the 


outcome of the spectral luminosity computations of our model 
(as dehned by eq. 0 ) for the case of NGC 5548 that we have 
chosen as a test of the model itself. 

The specihc meaning of the various curves in Fig.^is ex¬ 
plained hereafter. The soft photon (IR-Opt.-UV) spectral en¬ 
ergy distribution resulting from the observational data in the 
literature outlined above for NGC 5548 is shown in Fig. ^ as 
the continuous curve extending up to £ = 0.1 keV. The dotted 
curves represent the contribution to our model spectrum due to 
inverse Compton emission at three different times f,- = 0.02, 
1.0, 5.0 s (after the high energy electron injection in the loop), 
computed from expression 0 . It is apparent that the spectral 
energy distribution strongly depends on time and, in the X-ray 
energy range, it decreases with time, due to the evolution of 
the relativistic electron population, i.e. to its depletion starting 
from the high energy end of the distribution. 

For the results shown in Fig.^ the initial physical parame¬ 
ters of the relativistic electron distribution injected in the loop 
at f = 0.0 s are those specified in the caption of Fig. [Q Here 
we introduce as a significant parameter the product WoVioops, 
with Mo the normalization constant in relation ([0 defining the 
electron energy distribution, and Vioops the total volume oc¬ 
cupied by the active loops; its physical meaning is such that 
(noVioopsdj) represents the number of accelerated electrons in 
the whole flaring volume in an energy interval dy. In ad¬ 
dition, the other two physical parameters entering the emis¬ 
sion computations (see eq. ®) have been chosen as £uv-x = 
1 X lO''^ cm, Mth = 3 X 10^ cm“^; we note in passing that 
these values are those for which the best fit of the spectrum 
of NGC 5548 is obtained, as described in the following. 

There is one more curve appearing in Fig. [Q i.e. the con¬ 
tinuous one extending over the whole X-ray energy range, and 
in the following we discuss its meaning and the way it is com¬ 
puted within the framework of our model. From Fig. [J it is 
apparent that the computed emission changes on short time 
scales. Especially around a few hundred keV the decrease of 
the spectral energy distribution value with time is significant, 
owing to the depletion of high energy electrons due to inverse 
Compton and synchrotron losses. Note that we have dehned the 
quantity femis as the time the spectrum at £ = 10 keV takes to 
decrease by a factor of twenty with respect to its initial value. 
Eor the case shown in Fig.^a value of femis around 4.0 s can be 
derived. It is evident that every observational procedure which 
takes a time > femiss to get a spectrum will basically make a 
temporal mean of the time dependent spectra in the data acqui¬ 
sition. If X-ray emission in AGN is due to such a mechanism 
the fast evolving Hare spectra imply that a rapid succession of 
Hares must take place over different parts of the disk so that 
when one flare is fading out another one starts brightening. If 
the decay time for a flare is femis, we need one flare to start at 
least every femis seconds. We assume here a uniform distribu¬ 
tion in time of flares in order to reproduce the observed emis¬ 
sion. This would not be the case for a real AGN where the 
observed variability can be due to a non uniform distribution in 
time and in intensity of the flares. However, this scenario does 
not reduce the applicability of the model, since the observa¬ 
tional procedure does introduce an averaging operation on the 
emitted spectra on longer time scales. In conclusion, what we 
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observe is a temporal mean of the emitted spectra, i.e., in terms 
of spectral luminosity, 

iMEAN(£’) = - - I L\c(E,t)dt. (12) 

^emis Jo 

The result of this temporal mean in terms of an energy spec¬ 
trum (£’/(£’)mean) is drawn in Fig. [2 as the continuous curve 
extending up to high energies. 

The fast time evolution of the 1C emission has important 
consequences. In fact, since the estimated loop length scale is 
(see Sect. 2) i/flare > 3/?s ~ 6.2x10^^ cmforMBH = TxIO^Mq, 
the following relationship holds 

Gyn ~ ^flare/t^ ^ 2.07 X 10 S > temis- 

The numerical evaluation above refers to the specific case 
of NGC 5548, but, when repeating the same estimate for a 
lower value of the central black hole mass, say IO^Mq, we get 
the same ordering of the relevant timescales, namely fdyn - 
77flare/c > 300 s > femis- Therefore, we can generally con¬ 
clude that we do not need to take into account that at each 
electron reflection inside the magnetic loop a certain fraction 
of the electron beam precipitates into the disk, as explained at 
the beginning of Sect.2. 


4.3. Parameters of the model and their role 


The mean spectrum shown in Fig. [^depends on many free pa¬ 
rameters: the scale length of the corona, /?uv-x, the electron 
distribution parameters, yoi, 7 o 2 , WoFioops and s, the parameters 
entering the electron losses, namely the magnetic field B and 
the plasma thermal density nth- As for the other scale length ap¬ 
pearing in the equations relevant to the definition of the model 
spectrum, that is the characteristic linear size of each sin¬ 
gle emitting loop region, in our scenario it is directly related to 
the global corona scale length parameter defined above, namely 
^uv-x, by means of our specific choice of simplified geomet¬ 
rical structure; in fact, supposing that approximately a fraction 
w = 1 /10 of the total coronal volume is flaring, i.e. is responsi¬ 
ble for X-ray emission, at any given time, and that the number 
of simultaneously active loops (emitting flares) is q, the rela¬ 
tionship between the global coronal region dimension and the 
scale length of each emitting structure (flare) /?em is the follow¬ 
ing 


- Rvy-X 



(13) 


This choice is, of course, somewhat arbitrary, but nevertheless 
in its simplicity sufficiently general; in the fitting procedure de¬ 
scribed and discussed in the following Sections we have chosen 
for the number q of simultaneous flares at any given time a rep¬ 
resentative value of 10, as inferred by Haardt et al. J1994t . but 
this parameter can be changed and the related consequences are 
briefly discussed in Section 6 . 

The IR-UV spectrum which influences both the electron 
losses (see eq. 0 and 0) and the IC emission (eq. 0) is 
given as known for the time being, even though its time vari¬ 
ability makes the reliability of this assumption not obvious 


when simultaneous UV and X-ray observations are not avail¬ 
able. We have computed the mean spectrum, like the one shown 
in Fig. n changing in turn one of the above listed parameters. 
The effects each of them has on the mean spectrum are the fol¬ 
lowing: 


a) changing Ruv-x changes the shape of the mean spec¬ 
trum in a complex way, since this parameter appears not only 
explicitly in eq. 0 , but also, implicitly, through the evolution 
of the electron energy distribution, which depends on 1C losses 
(see eqs. 0 , 0 ); hence, the value of Ruv-x is determined by 
the comparison with the observations within a suitably selected 
range; 

b) a decrease in yoi makes the spectrum peak (which is 
around 100 keV in Fig.^ computed for yoi = 50) shift towards 
lower energies; 

c) changing 702 does not have any significant effect on the 
temporal mean of the spectrum for E < 300 keV, since the 
high energy portion of the electron distribution is in any case 
depleted very fast; 

d) MoFioops is essentially the normalization parameter of the 
resulting spectrum for £ > 0.1 keV (in other words, it is some¬ 
what proportional to the energy spectrum level) and the com¬ 
parison with observations determines its value; 

e) an increase in the value of s makes the spectrum steeper 
for £ > 100 keV; 

f) B has a complex influence on the mean spectrum; how¬ 
ever equally good fits of the data can be recovered for different 
values of the magnetic held intensity; 

g) nth does not change the shape of the mean spectrum for 
£<10^ keV. 


In conclusion, the parameters which determine the spec¬ 
trum in the range of BeppoSAX observations (0.1-200 keV) 
are £uv-x, «oFioops, Toi and s. In our study these parameter 
values will be defined by the fit to the observational data of 
BeppoSAX for NGC 5548. On the other hand, the values of 
the other three parameters 702 , B and n± are chosen by us be¬ 
fore the fitting procedure. We set 702 = 1000, since we have 
verified that this has no influence on the physics of the model. 

As for the other two input parameters, B and nth, the fol¬ 
lowing general considerations help us to estimate reasonable 
values. Assuming a value ^ 7 x 10^ Mq for the central black 
hole mass of NGC 5548 (Wandel et al. I2000> . the resulting 
Schwarzschild radius is Rs ^ 2.07 x 10'^ cm and an estimate 
of the lower limit for the length scale of coronal structures can 
be derived as £uv-x > 3 Rs - (£uv-x)inm (see Sect. 2). This 
lower limit for the length scale of the corona implies in turn 
an upper limit for the thermal density of the background ther¬ 
mal coronal plasma. In fact, as already mentioned in Sect. 2.1, 
it is our requirement that the thermal plasma in the coronal 
structure does not affect significantly the resulting high energy 
spectrum; this can be translated into a condition on the optical 
depth to scattering of the thermal component, tj, that must be 
tt = WthO^T^uv-x << 1. To meet this condition a limitation on 


1 


I 


the thermal density is obtained: n± « ^ 

2.4 X 10^° cm“^. We just note here 


-1 


0'T(/?UV-X)m 


(crT6.2 X 10^^ cm)- 
that for lower Mbh we expect a scaling of the limit on n±, but 
for Mbh ~ lO^M© it would anyway be nth 10** cm^^. In 
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our model we choose Wth = 3x10* cm~* which substantially 
fulfills the condition above, tt « 1, on a linear length scale 
^uv-x < 10'^ cm and, as we shall see in the next Section, 
“a posteriori”, certainly is in accordance with that same condi¬ 
tion for the value of R\jv-x - 1 x lO^'* cm that we find from 
the fitting procedure on NGC 5548. It is important to mention 
here that no specific evaluation of the temperature of this ther¬ 
mal component has been done, nor it is the aim of this work 
to do it. We only can say that, since we assume that most of 
the energy discharged in the coronal loop goes into accelerat¬ 
ing the relativistic electron population, we do not expect (ac¬ 
tually we can exclude) that this temperature can be as high 
as the values estimated in the thermal Comptonization models 
for the X-ray coronal emission. On the contrary, it is plausi¬ 
ble that the thermal plasma in the loops thermalizes at com¬ 
paratively low temperatures, with respect to those typical of 
thermal Comptonization models for X-ray emission. This in¬ 
ference, even though purely qualitative, together with the re¬ 
quirement expressed above of small optical depth to electron 
scattering of the thermal plasma, is the reason why we can ne¬ 
glect any thermal component effects [i.e., no thermal comp¬ 
tonization effects, no recoil distortion of the high energy end of 
the X-ray spectrum (see Krolik lT999H . 

As far as the magnetic field is concerned, its value is limited 
by energetic considerations. In fact, in the working hypothesis 
that magnetic energy conversion is the source of the accelerated 
electron energy, the following relationship must be fulfilled 

g2 njoi 

3 — > Mo ToCro - 1 ) 'djo- 

In this expression the electron distribution parameters are de¬ 
rived from the fit, as described above, and hence an iterative 
procedure is required to find a suitable value for the magnetic 
field, since B itself enters the fit. 

For NGC 5548 we have found that values of B in the range 
~ 110 - 500 gauss can produce model spectra that do fit the ob¬ 
served X-ray spectrum, provided that the typical coronal scale 
length /?uv-x is suitably chosen: it turns out that smaller coro¬ 
nal structures must be associated with stronger magnetic fields, 
in order to produce similarly fitting model spectra. 


5. Comparison with observations of NGC 5548 

The very extended energy range of BeppoS AX data makes 
it extremely well suited for comparisons with our model. 

We have thus computed the mean luminosity spectrum, 
Lmean{E), for the case of NGC 5548 using the IR-UV param¬ 
eters listed in the previous section. Then we have generated a 
grid of our spectra using the XSPEC atable format in order to 
compare our model to observations. 

NGC 5548 was observed by BeppoS AX from 2001 July 
08, 08:35 UT, to July 11, 00:16 UT. We used BeppoSAX 
data obtained with the Low-Energy Concentrator Spectrometer 
(LECS; 0.1-10 keV; Parmar et al. li^Tt . the Medium-Energy 
Concentrator Spectrometer (MECS; 1.65-10 keV; Boella 
et al. ll997bT l and the Phoswich Detection System (PDS; 15- 
300 keV; Erontera et al. ll997t . The net exposures in the LECS, 


MECS and PDS instruments are 42 ks, 98 ks and 48 ks, re¬ 
spectively. The net count rates for the LECS, MECS and PDS 
instruments are 0.37 cts s“', 0.55 cts s“* and 0.90 cts s“* re¬ 
spectively. 

The BeppoSAX data were reduced using the SAXDAS 
2.3.0 data analysis package. Good data were selected from in¬ 
tervals when the elevation angle above the Earth’s limb was 
> 4° and when the instrument configurations were nominal. 
The standard PDS collimator dwell time of 96 s for each on- 
and off-source position was used together with a rocking an¬ 
gle of 210 arcminutes. LECS and MECS data were extracted 
using radii of 8 arcminutes and 4 arcminutes, respectively. 
Background subtraction for the LECS and MECS were per¬ 
formed using the standard background files. The PDS back¬ 
ground was estimated from the offset field according to the 
standard procedure. 

The LECS and MECS spectra were rebinned to oversample 
the full width half maximum (EWHM) of the energy resolution 
by a factor 3 and to have, additionally, a minimum of 20 counts 
per bin to allow use of the statistic. The PDS data were re¬ 
binned using the standard logarithmic binning recommended 
for this instrument. Data were selected in the energy ranges 
0.13-2.4keV (LECS), L65-10keV (MECS), and 15-200keV 
(PDS), where the instrument responses are well determined and 
sufficient counts obtained. 

5.1. Fitting procedure 

First of all we have fitted these data with our model - 1 - Fe Ka 
line-H Compton reflection hump in the energy range 3-^200 keV, 
by allowing the following free parameters to vary: Ruv-x and 
HoPioops- The fitting process has been repeated for different val¬ 
ues of the parameters yoi, s, in order to determine for which 
values the fit on Ruv-x and WoFioops best reproduces the spec¬ 
trum (see Section 5.2). The reflection hump component is mod¬ 
eled by appropriately using pexrav in XSPEC (where the re¬ 
quired input power-law with cutoff parameters are derived by 
recursively “fitting” this simple representation of the X-ray pri¬ 
mary spectrum to our IC model spectrum). The narrow Fe K 
line is simply modeled with a zgauss additive component. 

We start the fitting procedure over this limited energy 
range, excluding the “low” energy portion of the data, because 
we want at first to disregard the possible effects of a warm ab¬ 
sorber, which, at low energies, could easily hide our model’s 
spectral features. Over this energy range we obtain a very good 
fit (xl = 0.78 (d.o.f. = 88)), once we have chosen the values 
of the parameters of the injected electron distribution at time 
t=0, i.e. 7oi = 50, 7 o 2 = 1000, s - 3 (see Sect. 5.2). Fig. |2l 
shows a comparison of our best fit model to the data in this 
“medium-to-high” energy range. For the physical parameters 
directly defining the generation of our non-thermal IC model 
spectrum, we obtain the values listed in Table 1; the results of 
the fit for those parameters referring to the other components 
of the observed X-ray spectrum, such as the Iron line and the 
reflection hump, are not included in Table 1, and are described 
later. The parameter Wot appearing in Table 1 is defined as 

£ 02 

(y — 1) dy ; 
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Fig. 2. BeppoS AX spectrum of NGC 5548 and the best fit with 
our model + Fe Ka + reflection hump, in the range 3-200 keV 
(see text). 


this parameter represents the total number of relativistic elec¬ 
trons in the coronal structure volume and it is a function of 
{«oFioops, 701 , 702 , ■s)- We want to stress again that the^^ value 
in this case is very good, confirming that our model gives a 
good interpretation of the data in this energy range, which ba¬ 
sically corresponds to the energy interval over which the intrin¬ 
sic IC model strongly influences and determines the overall re¬ 
sulting spectral shape, together with the reflection component, 
which, in turn, depends on the intrinsic spectrum itself. 

As already noted above. Table 1 only shows the parameters 
of our IC model spectrum, but of course the fit determines also 
the parameters referring to the other components of the result¬ 
ing X-ray spectrum that we included, that is the Fe K line and 
the reflection hump. Our best fit selects a central energy for the 
line = 6.37 + 0.08 keV and a line width cr < 170 eV. As for 
the reflection hump, the resulting reflection factor R = Q/27r, 
where Q is the solid angle subtended by the reflecting matter 
(see Magdziarz & Zdziarski ll995t . is 0.8, and the cosine of the 
inclination angle of the reflecting slab-like material with the 
line of sight has been fixed at 0.9, like in Pounds et al. (I2003t . 
and Perola et al. (I2002t . It is worth noting that, if the cosine of 
the inclination angle is left free to vary, a value = 0.80 (87 
d.o.f.) is obtained, corresponding to a value 0.8 for the cosine 
itself and again R = 0.8, with the parameters reported in Table 1 
essentially unchanged. 

Table 1. The 3-200 keV physical parameters for our model, 
selected as described in the text. 


7oi 

= 50 

7o2 

= 1000 

5 

= 3.0 

B 

= 500 gauss 

Ruv-X 

= (1.0+"“) X 10*^ cm 

Mo, 

= (5.34);°g2) X 10“*’ electrons 

xl 

= 0.78 (88 d.o.f.) 



Fig. 3. Comparison of BeppoS AX data over the whole energy 
range observed with the IC “primary” spectrum + Fe line -H 
reflection component as derived from the 3-200 keV best fit; 
the necessity of some absorption component is apparent from 
the plot; see the text for further explanation. 


We then use these same parameter values, by fixing them 
in the XSPEC procedure, to compare the model (that we have 
found well fitting the 3-200 keV range) to the whole energy 
range of BeppoS AX data; the result of this comparison is much 
worse, giving an unacceptably large value of xl (~ 16). In 
fact, keeping fixed the parameters selected in the restricted en¬ 
ergy range (3-200 keV) and “freezing” the normalization con¬ 
stant of the LECS data group to a value such that the ratio of 
LECS/MECS normalization lies within the acceptable interval 
(see Eiore et al. rr^ and, at the same time, such as to make 
the model match the data in the low energy range (E $ 0.4), 
the model flux f(E) (keV cm“^ s“' keV^') that we obtain is in¬ 
stead significantly higher than the unfolded observed spectrum 
in the range 0.4- 3 keV; this is shown in Eig.|3 

This energy range is exactly the one where the expected 
effects of a warm absorber are most evident (Nicastro et al. 
I1999a> . Indeed, recent observations of this energy range with 
the unprecedented spectral resolution of XMM-Newton, have 
shown an astounding complexity of spectral features, emission 
and absorption lines, possibly blueshifted, (see Steenbrugge et 
al. 120031 Kaastra et al. I2002> revealing how complicated and 
composite a realistic description of the absorbing gas must be, 
certainly falling beyond the scope of the present work. In the 
context of Beppo-SAX data analysis in the “low-energy” spec¬ 
tral region, characterized by a limited spectral resolution (with 
respect to that of the grating spectrometers on board of XMM), 
similarly to what is done by Pounds et al. J2003t . we choose to 
give an approximate description of the warm absorber effects 
with a series of absorption edges, namely 5 edges. 

We postpone the discussion of these edges until Section 5.3 
and at present we only say that including them in the fitting 
procedure (each of them represented by a multiplicative com¬ 
ponent ZEDGE in XSPEC), and keeping the IC model spectrum 
parameter values fixed as those in Table 1, now we obtain a 
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Fig. 4. BeppoS AX spectrum of NGC 5548 and the best fit with 
our model + Fe Kor + reflection hump + absorption edges (see 
text). No soft excess is included. 


much better ht (x'y - 0.96 (133 d.o.f.)) over the entire observed 
energy range 0.12-200 keV; this is shown in Fig.|3 

5.2. Physical parameters for the generation of the 1C 
spectrum 

Table 1 shows the best fit values for the parameters defining the 
physical scenario that we have devised for the generation of the 
primary X-ray spectrum. 

We have also tried to fit models corresponding to different 
values of s and yd; the ones we have hnally chosen, given in 
Table 1, are those for which the best fit is obtained, i.e. corre¬ 
sponding to one of the lowest obtainable values. However, 
it is interesting to note that variations of these two parameters 
within ~ 15% of the reported values result in a ht which is 
still a reasonable one, if not slightly better; for example, for 
i = 2.5 the hnal^^ value turns out to be even lower than the 
one found for s - 3 (see Table 1). Nevertheless, we have pre¬ 
ferred and shown the ht obtained for the model corresponding 
to s = 3, since in this case the decrease of the spectrum at high 
energies is steeper than in the case corresponding to s = 2.5. 
Regarding the energy range of the injected electron distribu¬ 
tion, recent work on electron acceleration due to magnetic re¬ 
connection related processes (Lesch & Birk I 19971 Schopper et 
al. lD^ in the AGN context and, even more specihcally, in the 
AGN corona context, has shown that such processes guarantee 
acceleration up to maximum energies around y ~ 2000 , which 
is essentially in accordance with our choice of the upper limit 
of the energy range of the initial injected electron distribution. 

The values of Ruv-x and Atot reported in Table 1 have 
been directly derived from the best ht procedure. In the hy¬ 
pothesis that approximately one tenth of the coronal volume 
is Haring at the same time and that the emission is due to a 
series of q - 10 almost simultaneous Hares, the characteris¬ 
tic linear scale of a typical emitting region can be estimated 
as Rem ^ Ruv-x/(10<7)'/3 ^ 2.15 x lO^^ cm and the corre¬ 

sponding minimum variability time scale is ~ 7.2 x 10^ s, in 
accordance with the the data reported by Chiang et al. (EOOOj. 


In this framework the density of non-thermal electrons can be 
evaluated as 


Wrel 


Mot 

10(47r/3)RL 


^ 1.28 X 10^ cm-^ 


and the associated optical depth to (non-thermal relativistic) 
electron scattering is Tjei = ni-eiCTRem - 1.8 x 10“^. With such 
a low optical depth, secondary inverse Compton scattering ef¬ 
fects can be safely neglected. The computation of how many 
double scattered photons attain E - 300 keV shows that their 
number is signiHcantly lower than that of the photons that have 
been scattered only once. 

Finally, we want mention that a similarly good Ht from our 
model can be obtained also for larger values of Ruv-x, and, 
correspondingly smaller magnetic Held values, as noted at the 
end of Section 4.3. In particular, for B - 120 gauss, from the 
Ht in the 3-200 keV range, we obtain = 0.76 (87 d.o.f.), for 
Ruv-x = (6.01°;“)x 10'4 cm and Mot = (1.40+°;°|) x 10^' elec¬ 
trons, with essentially the same resulting parameters for the 
iron line, and with R = 0.8 and (0.67^g | 3 ) for the cosine of 
the inclination angle for the pexrav component, deHning the 
reflection component of the spectrum. Extending the Ht to the 
whole energy range of BeppoSAX data in the same way de¬ 
scribed in the previous section, we obtain xl - 0.95 (132 
d.o.f.), thus giving as good a Ht as the one described in the 
previous Section. For this case Rem - 1.29 x 10'^ cm, so that 
the relativistic electron density is Mrei - 1.56 x 10^ cm“^, and 
the associated optical depth to scattering is Trei - 1.3 x 10^^. 
Nevertheless, we have chosen to fully illustrate the Ht corre¬ 
sponding to Ruv-x = 1 X lO'"^ cm in the previous and in the 
present Sections, following the outcomes of the discussion of 
Section 2 regarding the estimate of the typical scale length of 
active loops, which indicate that smaller size coronal loops are 
favoured. 


5.3. Other components of the observed X-ray 
spectrum 

We want to mention here that the Ht we have obtained seems 
to indicate that no signiflcant soft excess is required from the 
2001 data we have analyzed. This might seem in disagreement 
with Pounds et al. (I2003t (who detect an “apparently unam¬ 
biguous excess of soft X-ray flux below 0.7 keV”); however, 
this is only apparent, since their conclusion comes from the use 
of a simple power-law to describe the primary spectrum. On the 
contrary, in our case the X-ray primary spectrum itself is the 
result of inverse Compton reprocessing of soft photons by our 
evolving non-thermal distribution of relativistic electrons, and, 
as such, turns out to be characterized by a more complex be¬ 
haviour. Speciflcally, it shows a slight hump right in the energy 
region around and below ~ 1 keV, as it is suggested by Fig. ^ 
and could be better seen with a more appropriate scale. This 
is just the time-averaged effect of the evolution of the lower 
limit of the electron energy distribution, ji = 71 (f) < 7 oi, and 
it turns out to be essentially enough to suppress the necessity 
of an additional soft excess component. We note, however, that 
we might indeed “force” a soft excess component in the frame¬ 
work of our model and still obtain a reasonable Ht; in other 
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words, we do not exclude the possibility of a soft excess com¬ 
ponent, but we stress again that it turns out to be not necessary 
to introduce it within our modeling scenario. 

Let us briefly discuss the XSPEC components that we have 
added to our resulting IC spectrum. In fact, we have introduced 
these components since it is well known that they do have an 
influence on the “low” X-ray energy range for what regards the 
warm absorber, and in the medium-to-high energy range for 
the Iron emission line and the reflection hump (Nicastro et al. 
120001 Petrucci et al. I2000> : as a consequence, to fit the whole 
BeppoSAX energy range it is necessary to account for them. 
However, we want to stress that we are not interested here in 
discussing the nature/origin of these reprocessing components 
in detail, but instead we want to put emphasis on the feasibility 
of an “intrinsic” emission model such as the one we have envis¬ 
aged. Furthermore, the parameters we obtain from the general 
fit for the warm absorber and for the other reprocessing com¬ 
ponents are not in disagreement with the values cited in pre- 
XMM analyses in the literature (Nicastro et al. 120001 Pounds 
et al. l200^ . 

As for the five absorption edges through which we give 
a simplified description of the warm absorber effects, four of 
these edges have been fixed to the theoretical OVII, OVIII and 
NeIX/X and MgX/XI edge energy values (and these are essen¬ 
tially the four edges mentioned by Pounds et al. J2003H . The 
fifth edge, which is required if we want to obtain a good fit, 
turns out to be placed at ^ 0.5 keV, and it could correspond 
to the deep absorption edge of CIV predicted by photoioniza¬ 
tion models of the warm absorber (see Nicastro et al. I1999al 
I1999bt . However, what is of interest to our present purpose is 
that we do obtain a good fit by allowing for absorption edge 
energies and resulting optical depth values which are in accor¬ 
dance with published studies of warm absorbers prior to the 
high energy resolution observations of XMM. In fact, this al¬ 
lows us to devote our attention to the “primary” X-ray contin¬ 
uum as emitted in the non-thermal scenario we have devised. 
Moreover, it is not very meaningful any longer to dwell on the 
meaning and origin of “low-energy” spectral features as they 
appear from BeppoSAX data, since CHANDRA and XMM, al¬ 
lowing for much higher spectral resolution and therefore much 
more detail, have now provided a much more complex picture, 
revealing a richness of lines and spectral features whose anal¬ 
ysis requires a different approach and different means of study, 
and, possibly, showing no evidence at all for the oxygen edges 
previously identified in lower energy resolution observations 
(Kaastra et al. l2000t . 

6. Application of the model to AGN X-ray sources: 
general considerations 

In the previous section we have shown that the X-ray spec¬ 
trum of the specific Seyfert 1 source NGC 5548 can be interest¬ 
ingly well explained by the non-thermal and intrinsically non¬ 
stationary model that we have devised. Now we turn our atten¬ 
tion to the analysis of the descriptive capability of the model 
with respect to a few general observed trends and properties 
of Seyfert I’s X-ray spectra. We do not get into the detailed 
and quantitative description of a fit, like we did in the previ¬ 


ous two sections; in fact, since in our model the reprocessed 
spectrum depends on the input soft radiation spectrum, and 
on the loop magnetic field through the locally generated syn¬ 
chrotron component, which represents a seed component for 
inverse Compton reprocessing as well, to compare its outcome 
with the average Seyfert 1 X-ray spectrum, we should have de¬ 
fined an averaged opt.-UV spectrum to be reprocessed and an 
average magnetic field, and then fit the resulting X-ray spec¬ 
tral distribution to the “averaged” observed spectrum, obtaining 
values of the physical parameters coming into play that would 
define “average” physical properties of the emitting regions. 
As a matter of fact, especially because of the significant depen¬ 
dence of these parameters on the physical condition of the very 
central region of the AGN, i.e. on the effective value of the cen¬ 
tral black hole mass ultimately, a procedure of this kind would, 
in our opinion, lead to not very meaningful results. On the con¬ 
trary, we point our attention to the capability of our model to 
reproduce a very general spectral characteristic of the observed 
Seyfert X-ray spectra, that is the relatively small dispersion of 
the distribution of the spectral index a of a power-law like de¬ 
scription of the “primary” spectral flux (/(£) oc at least in 
the energy range ~ 2 -10 keV around an average value a ~ 0.9 
(Nandra & Pounds 179941 Nandra et al. ll997al Zdziarski et al. 
119^ Zdziarski et al. 120001 Gondek et al. [19961 Matt EOOTI for 
BeppoSAX broad-band observations), by now considered as a 
“canonical” value. It is appropriate to note here that, strictly 
speaking, such an inference on the primary spectrum is some¬ 
what model dependent, since the “canonical” value a - 0.9 is 
obtained fitting the observed spectrum with a series of compo¬ 
nents (primary and reprocessed, see also Section 5.1) that are 
a priori supposed to be descriptive of the general mechanisms 
at work in the definition of the resulting X-ray radiation. In 
particular, the existence of a reflection component turns out to 
be determinant; for an interesting, and in a way “historical”, 
analysis of the developments of the global phenomenological 
understanding of X-ray spectra of Seyfert 1 nuclei, we refer 
for example to Nandra & Pounds’s J1994t work on GINGA 
data, showing that the strong flattening of the spectrum in the 
10-18 keV range, with respect to the 2-10 keV behaviour, was 
well accounted for by introducing a reflection component in 
the spectral fitting, thus giving a relatively narrow distribution 
of values of the spectral index of the power-law description of 
the primary spectrum around an average value ~ 0.9. 

Back to the point, it is our purpose to show that our model can 
quite generally reproduce a spectrum that fairly well approx¬ 
imates a power-law behaviour in the range 2-10 keV with a 
spectral index that does match the relatively narrow range of 
values around a - 0.9 inferred from the general interpretation 
of Seyfert X-ray spectra observations. In fact, we restrict our 
analysis to the energy interval delimited by ~ 2 and 10 keV re¬ 
spectively because of two orders of reasons, which we explain 
in the following and are anyway intimately connected with 
each other. Indeed, on one hand our model primary spectrum, 
when considered on a much broader X-ray energy range (such 
as the one observed by BeppoSAX), necessarily (i.e. due to its 
origin) is much more complex than a simple power-law. On the 
other hand, within the general interpretative framework, with 
the exception of the generally narrow (see Bianchi et al. l2004t 
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Iron line at ~ 6.4 keV, the effects of the reprocessing compo¬ 
nents usually accounted for in the theoretical explanation of the 
resulting broad-band X-ray spectrum are expected to be much 
more conspicuous outside of this ~ 2 - 10 keV range, (because 
of absorption on the lower-energy side, and of reflection on the 
higher-energy side). Therefore, outside of the 2 -10 keV range, 
deviations of our model spectrum from a simplistic power-law 
behaviour are mixed up with reprocessing features/distortions 
and can be disentangled only by altogether fitting the observed 
spectrum, possibly giving sort of different results for the re¬ 
processing component parameters (with respect to the simple 
primary broad-band power-law description) but in no way re- 
jectable on a theoretical basis. Keeping in mind the discussion 
above, to attain our goal, we have chosen the following proce¬ 
dure. First, we had to define some sort of “average” properties 
of a typical Seyfert 1 spectrum in the optical-UV range. Several 
works in the literature attempt to define these properties for the 
general category of radio-quiet AGNs, or for QSOs only (see 
for example Risaliti & Elvis 120041 and references therein), but, 
in order to select an average opt.-UV spectral behaviour more 
specifically referring to the class of Seyfert Is, we referred to 
Koratkar & Blaes’s work, specifically to their Fig. 1, in 

which the authors show an average SED (spectral energy dis¬ 
tribution) for Seyfert Is, distinguishing it both from radio-quiet 
quasars (QSOs) and radio loud ones. Approximately fitting a 
opt.-UV Sy-1 spectrum from Koratkar & Blaes’s J1999t Fig. 1, 
we computed a grid of resulting primary X-ray spectra from 
our model, changing the physical parameters that we found 
more significant in the detailed fitting of NGC 5548’s X-ray 
BeppoSAX spectrum as performed in Section 5, namely Ruv-x 
and Atot, together with the magnetic field in the active region, 
B. We find that, for a reasonable range of Atot values, and for 
a given Ruv-x, we can quite easily define a range of values of 
the magnetic field B corresponding to which our model spec¬ 
trum in the range 2-10 keV can be well described as power-law 
like, with spectral index a values between 0.85 and 0.95. To 
quantify our assertion, for example, this typical slope range for 
the power-law best describing our spectrum can be obtained 
for B ^ 100 - 200 gauss, when Ruv-x = 6 x lO^'* cm, whereas 
for Ruv-x = 1 X lO^'* cm this same a-range is reproduced for 
B ^ 400 - 600 gauss. This is in accordance with the trend 
that we have briefly discussed at the end of Section 4.3, refer¬ 
ring to the specific fitting of NGC 5548’s BeppoSAX spectrum; 
in fact, we do find that, even “starting” from an an average 
opt.-UV spectrum to be IC reprocessed by the coronal rela¬ 
tivistic electrons, an increase of the coronal structure extension 
requires smaller values of the magnetic field to reproduce the 
canonical power-law slope, i.e. to better fit observed spectra. 

Another issue we have to confront with is X-ray variability, 
its general behaviour and its relationship with variability prop¬ 
erties at other wavelengths. X-ray spectral variability studies 
have shown that many Seyfert spectra soften when the contin¬ 
uum flux increases (Nandra et al. ll997bl Markowitz & Edelson 
120011 Markowitz et al. l2003t . Also in NGC 5548 temporal vari¬ 
ability has been observed at several wavelengths; in particular, 
UV-EUV and simultaneous X-ray variability data are available 
(Chiang et al. 120001 Haba et al. 120031 Kaastra et al. l2004t . and 
the same trend, namely an UV-EUV intensity increase larger 


than the X-ray one, is present. Another characteristic behaviour 
which appears both in NGC 5548 (Chiang et al. 120001 Kaastra 
et al. l2004> and also in some other sources (see, e.g., McHardy 
et al. 120041 Blustin et al. I2003II is the fact that hard X-ray in¬ 
tensity variations show a delay with respect to the soft X-ray 
ones. 

A detailed comparison of our model results with the ob¬ 
served variability properties is out of the purpose of this pa¬ 
per. However, it is important to understand whether our model 
could reproduce the most relevant variability trends outlined 
above. In our model both an increase in the magnetic field in¬ 
tensity and a different geometry, i.e. less numerous but more 
extended active loops, can easily reproduce a softening of the 
spectrum. More specifically, an increase of the magnetic field 
intensity of 25% can induce an increase in the observed spec¬ 
tral luminosity of 22% at 0.3 keV and of 7% at 3 keV, while a 
decrease of the number of active regions from ten to three can 
induce an increase in the observed spectral luminosity of 28% 
and 10 %, respectively corresponding to the same two energy 
values cited above. A non-homogeneous magnetic field show¬ 
ing different geometrical configuration and intensity in differ¬ 
ent regions of the accretion disk, owing to the disk rotation 
and to its intrinsic temporal variability, would produce changes 
such as the ones mentioned above in the coronal magnetic 
structure and in the number and properties of flaring loops. 
This picture is in accordance with the idea of Kaastra et al. ( 
I2004t regarding the possible presence of a large hot active re¬ 
gion above the disk. In this context, magnetic field properties 
can be imagined to evolve in time in such a way to justify the 
observations of the hard X-ray part of the spectrum still in¬ 
creasing, while the soft one decreases, in accordance with the 
reported time lag of hard X-ray radiation. 


7. Discussion and Conciusions 

In this paper we have developed the idea that AGN coronae 
resemble the stellar ones, where flares powered by magnetic 
reconnection inject accelerated electrons into magnetic, loop¬ 
like, structures. In this framework, X-ray emission from AGNs 
is due to the non-thermal, time-evolving distribution of rela¬ 
tivistic electrons and the dominant emission mechanism is in¬ 
verse Compton scattering of IR-UV soft radiation emitted by 
the underlying accretion disk. The model has been tested on 
one of the brightest Seyfert 1 galaxies, NGC 5548. 

Despite making exceedingly simplifying assumptions, our 
model gives a very good fit to the X-ray spectrum of NGC 5548 
(see Figs. 121 and 13 as observed by BeppoSAX in July 2001, 
with physically reasonable values of the model parameters. In 
addition, the spectral index of the initial electron energy distri¬ 
bution, s, that we derive turns out to be quite in accordance with 
the values inferred for solar flares, hence supporting the idea 
that magnetic reconnection is related to the electron accelera¬ 
tion mechanism. Moreover, ourRuv-x - 4.8Rs is substantially 
consistent with the characteristic size of a coronal structure for 
an object like NGC 5548. If ten flares are active at the same 
time in the corona (see Sect. 5), the estimated size of the blob¬ 
like structure implies that significant variations in the observed 
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luminosity are on time-scales longer than ~ 7.2 x 10^ s, an 
estimate which is in accordance with X-ray luminosity varia¬ 
tions observed up to now (Chiang et al. 120001 Nicastro et al. 
mm . Since our model is based on a time-evolving emission 
process, possible variations of the number of active flares and 
of their intensity allow to decouple X-ray intensity variations 
from the UV ones and also to reproduce the observed changes 
in the emission spectrum (see section 6). Another advantage of 
our model is that it does not require the introduction of a soft 
X-ray emission component, since the model itself reproduces 
the so called “soft excess” (see Sect. 5.3). 

A qualitative comparison with general trends inferred from 
observations of the X-ray spectra of Seyfert 1 active nuclei has 
been discussed in Section 6: both the “average” spectral index 
a ^ 0.9 generally determined for power-law approximations of 
the “primary” spectrum and the variability patterns observed 
for several sources (and in particular for the one, NGC 5548, 
we analyzed in more detail) briefly described in the previous 
Section can be rather easily reproduced in the framework of 
our model. This is encouraging in what concerns the feasibility 
of the scenario we propose in the present work. 

As outlined in the introduction, two main different fami¬ 
lies of models have been developed in order to understand the 
physics of the X-ray emitting plasma; those in which the energy 
distribution of the electrons responsible for X-ray emission is 
thermal and those in which it is instead non-thermal. It is inter¬ 
esting to highlight the most important points distinguishing the 
presently proposed scenario from the above cited families of 
models. Differences with thermal models are straightforwardly 
identihed, since our model is based on the impulsive injection 
of highly relativistic electrons non-thermally distributed in en¬ 
ergy (and evolving in time due to their energy losses), interact¬ 
ing with seed soft photons to produce a time dependent X-ray 
spectrum through a single-scattering inverse Compton process. 
On the contrary, the usual thermal models are characterized 
by a steady-state trans-relativistic thermal energy distribution 
for the electrons, which interact with the soft radiation held 
again through the inverse Compton mechanism, but in a multi¬ 
ple scattering process, thermal Comptonization, with each one 
of the scatterings involving a relatively small amount of en¬ 
ergy exchange. In our scenario we do allow for an “underly¬ 
ing” thermalized coronal plasma component, but its physical 
properties differ from those of “thermal models” quite signih- 
cantly, as discussed in Sect. 4.3; in particular, the density of the 
thermalized component is much lower than that estimated for 
thermal models. As a consequence, since the optical depth to 
scattering of this component is much smaller than one, no sig- 
nihcant effect from the Compton interaction between thermal 
matter and photons can be expected, not even Compton “re¬ 
coil” (see Sect. 4.3). The “primary” spectrum is therefore only 
determined by the non-thermal relativistic electron interactions 
with the radiation held. 

To this respect, it is also important to stress the properties of 
our coronal emitting component, i.e., of the non-thermal rela¬ 
tivistic electron population. In fact, also this component is opti¬ 
cally thin to scattering, since its Thomson optical depth is very 
low (~ 0.0018 for the case of NGC 5548), and this justihes our 
choice of a simple single-scattering description. 


Concerning the “standard” stationary non-thermal models 
(see Svensson 119941 Ghisellini I1994t . two main differences 
with respect to our present model can be identihed and are dis¬ 
cussed in the following; however, it is important to notice that 
these same two differences also apply to the comparison with 
stationary thermal Comptonization models mentioned above. 

The hrst signihcant difference is the fact that our model 
is intrinsically non-stationary and the evolution in time of the 
resulting spectra (as it has been shown in Section 4.2) turns 
out to be a decisive factor in order to reproduce the observed 
steepening of the X-ray spectrum above 100 keV, in accordance 
with OSSE and BeppoSAX results. 

Secondly, the compactness parameter (see below for its ex¬ 
plicit dehnition, as given by Guilbert et al. I1983t . consistently 
evaluated for the coronal emitting regions in our scenario, is 
lower than the critical limit for pair production to be signihcant, 
thus avoiding, in our case, the problem of the annihilation line, 
generally expected from standard, stationary non-thermal mod¬ 
els, but not detected by OSSE (Johnson et al. 1 19971 Zdziarski 
et al. ll997t . In fact, following Guilbert et al. (I1983t . the com¬ 
pactness parameter of the coronal structure can be estimated as 
/corona ~ cTj!{msC^){LxlR\}\i-x), where Lx is the source X-ray 
luminosity and J?uv-x is the global extent of the region where 
the IC reprocessing occurs. Using for Lx values in the range 
2-3x10"^^ erg/s (as reported by Nicastro et al. 120001 Perola 
et al. 120021 Bianchi et al. 120041 for NGC 5548), we get an es¬ 
timate /corona ~ 4.9 - 7.4, reasonably below the well known 
limit for pair production to become effective, namely I ^ 4;7r. 
Even lower values can be estimated for the compactness pa¬ 
rameter for the single emitting loop (the “direct” analogue of 
the “blob” of Haardt et al. [T994J, since in that case we have 
/loop ~ o-T:l(meC^)(Lxlq)IRem ~ 2.3 - 3.4, to be compared with 
the reference value ~ 30 of Haardt et al. Therefore, in 

our framework neglecting pair production/annihilation effects 
on the spectrum turns out to be fully justihed from the resulting 
conditions of the emitting structure, and this supports the con¬ 
sistency of our treatment. It is also noteworthy that this condi¬ 
tion differentiates our model from the thermal Comptonization 
ones as well, as it is clear from the comparison above with com¬ 
pactness properties of the thermal “patchy-corona” model of 
Haardt et al. lfT9^ . 

Still referring to the issue of a comparison of our presently 
proposed scenario with preexistent models, we have to mention 
that models including both a thermalized (Compton efficient) 
electron population and a non-thermal relativistic component 
have been also devised; these models have been developed to 
account for the fact that a purely thermal (Maxwellian) distri¬ 
bution for energetic electrons is quite difficult to achieve (see 
Coppi fT^^ . 

On one hand thermalization mechanisms, which are more 
efficient than Coulomb collision for energy exchange, have 
been explored, although, according to Coppi J1999t . “realis¬ 
tic calculations are in general quite difficult, especially near 
an accreting black hole, where we still have relatively poor 
knowledge of the exact physical conditions”. It is far beyond 
the scope of the present paper to analyze this topic in detail; 
we just mention one interesting model belonging to this hrst 
sub-class, namely the so-called “synchrotron boiler”, devised 
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by Ghisellini et al. (fT^ and further developed by Ghisellini 
et al. (I1998t . based on efficient self-absorption of synchrotron 
photons, emitted by electrons in the source, by different elec¬ 
trons of the same population, thus “rapidly” exchanging en¬ 
ergy and producing relaxation of the particle population itself 
to a resulting steady-state quasi-Maxwellian distribution, with 
a high-energy tail. It is important to notice that the main as¬ 
sumptions for this mechanism to be efficient, and thus for ther- 
malization to effectively occur (Ghisellini et al. ll998t . are that 
i) the mean energy, < y >, of the energetic electrons is just a 
few, and ii) the magnetic energy density is dominant with re¬ 
spect to the radiation energy density, together with the further 
condition that the soft photons to be Comptonized only arise 
from the reprocessing of about half of the hard radiation by 
cold (disk) matter close to the active (Comptonizing) region. 
None of these conditions is met in our scenario, in which no 
such straightforward and stringent relation between the exter¬ 
nal soft radiation energy density and the hard (i.e. X-ray) one 
is required, thus allowing our magnetic energy density to be 
somewhat lower than the soft seed radiation energy density, 
and in any case not dominant (even though magnetic energy 
is in our view the reservoir for the electron impulsive energiza¬ 
tion, and, therefore, ultimately, for the high energy emission); 
moreover, contrary to assumption i), mentioned above, for the 
thermalization mechanism to be efficient, the accelerated elec¬ 
trons we consider are characterized by much higher electron 
energy » 1). 

On the other hand, models intrinsically allowing for a 
“lower” energy portion of the electron distribution of a 
Maxwellian form and a steady injection of high energy elec¬ 
trons in a power-law or delta-function distribution have also 
been devised (Coppi I1999t . These are the so-called hybrid 
thermal/non-thermal models and the resulting steady state dis¬ 
tribution is in fact a Maxwellian plus a high energy non-thermal 
tail, whose relative importance depends basically on the ratio 
of the power supplied to the thermalized component (heating 
rate) and the power injected in non-thermal electrons (fth/^nth 
in terms of compactnesses) as well as on the non-thermal injec¬ 
tion spectrum (especially the characteristic y range). As a con¬ 
sequence, the resulting high energy spectrum depends on these 
parameters as well. It is well known (see Ghisellini et al. ll993t 
that, as long as the radiative process is multiple Compton scat¬ 
tering and the maximum particle energy is a few MeV, the re¬ 
sulting high energy photon spectrum is indistinguishable from 
the thermal expected one, and is basically independent of the 
details of the distribution of particles. When the ymax of the 
non-thermal electrons is instead quite high, the spectral prop¬ 
erties will depend on the ratio but, even when this ratio 

is » 1, the spectrum should be distinguishable from the ther¬ 
mal one at gamma-ray energies. These models are at present 
mostly applied to galactic black holes in their soft states, that 
cannot be fitted by the thermal models, since the correspond¬ 
ing spectra do show power-law like tails at least up to 1 MeV 
(although relatively weak), and, as a consequence, appear to re¬ 
quire a high energy non-thermal tail in the electron distribution 
contributing to the formation of the spectrum itself; on the con¬ 
trary, for X-ray spectra of normal broad-line Seyfert 1 active 
nuclei, thermal Comptonization models have been preferred, 


up to now. It is appropriate, here, to cite a significant statement 
in Coppi’s JI999t review on hybrid models, again enlightening 
a substantial difference with our presently proposed scenario; 
the spectra produced by Coppi’s code for hybrid models are 
“steady state, while the real spectra that are being fit are typi¬ 
cally time integrations over many flares...”. 

In the light of all the arguments above discussed, we can 
conclude that while the comparison of our model with observa¬ 
tions suggests that our model is plausible and can be considered 
a valid alternative to thermal Comptonization models, a deeper 
analysis is still necessary to fully understand the physics of the 
emitting medium in AGN coronae. 
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